Naval  Research  Laboratory 

Washington,  DC  20375-5320 


NRL/MR/6707— 99-8348 


Parameter  Space  for 
Collisionless  RF  Sheaths 


Wallace  M.  Manheimer 
Plasma  Physics  Division 


April  8,  1999 


Approved  for  public  release;  distribution  is  unlimited. 


l 


19990413  055 


) 


OTIC  QUALITY  INSPECTED  8 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson 
Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503. 

1  AGENCY  USE  ONLY  ( Leave  Blank) 

2.  REPORT  DATE 

3.  REPORT  TYPE  AND  DATES  COVERED 

April  8,  1999 

Interim 

4.  TITLE  AND  SUBTITLE 

5.  FUNDING  NUMBERS 

Parameter  Space  for  Collisionless  RF  Sheaths 

JO  #67-7641 -A9 

6.  AUTHOR(S) 

Wallace  M.  Manheimer 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

Naval  Research  Laboratory 
Washington,  DC  20375-5320 

NRL/MR/6707-99-8348 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Office  of  Naval  Research 

800  North  Quincy  Street 

Arlington,  VA  22217-5660 

10  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 

12a  DISTRIBUTION/AVAILABILITY  STATEMENT 

12b.  DISTRIBUTION  CODE 

Approved  for  public  release;  distribution  is  unlimited 

13.  ABSTRACT  ( Maximum  200  words) 

As  plasma  processing  reactors  approach  higher  ion  density,  the  sheath  models  which  neglect  the  rf  response  of  the  ions  become 
invalid.  This  work  shows  that  at  arbitrary  ion  density,  the  nature  of  the  collisionless  rf  ion  sheath  can  be  described  in  a  number  of 
different  regimes  of  parameter  space.  These  regimes  can  all  be  visualized  on  a  single  two  dimensional  plot  where  the  horizontal  axis 
is  the  ion  plasma  frequency  divided  by  the  frequency,  and  the  vertical  axis  is  the  electron  oscillating  velocity  divided  by  the  ion  sound 
speed. 

14.  SUBJECT  TERMS 

15.  NUMBER  OF  PAGES 

18 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 

OF  REPORT 

18.  SECURITY  CLASSIFICATION 

OF  THIS  PAGE 

19.  SECURITY  CLASSIFICATION 

OF  ABSTRACT 

20.  LIMITATION  OF  ABSTRACT 

UNCLASSIFIED 

UNCLASSIFIED 

UNCLASSIFIED 

UL 

NSN  7540-01-280-5500 

Standard  Form  298  (Rev  2-89) 
Prescribed  by  ANSI  Std  239-18 

298-102 


1 


PARAMETER  SPACE  FOR  COLLISIONLESS  RF  SHEATHS 


I.  Introduction 

Plasmas  are  now  used  extensively  in  material  processing.  In  many  instances,  the  plasma 
supplies  a  stream  of  energetic  ions  to  the  surface;  and  this  stream  of  ions,  along  with  the 
various  neutral  chemistry  processes,  gives  rise  to  the  appropriate  surface  modification. 

For  a  steady  state  plasma,  the  passive  dc  sheath  (i.e.,  assuming  the  wall  at  the  floating 
potential)  gives  an  ion  energy  impinging  on  the  surface  of  typically  a  few  times  the 
electron  temperature.  Since  the  electron  temperature  is  typically  a  volt  or  two,  the  ion 
energy  may  be  perhaps  5  to  10  volts.  For  many  processing  applications,  higher  ions 
energies  are  required.  These  are  often  generated  by  biasing  the  substrate  with  an  rf 
current.  (If  the  substrate  is  an  insulator,  an  rf  potential  is  the  only  way  to  bias  it  in  a 
controlled  manner.)  The  rf  bias  gives  rise  to  a  sheath  which  accelerates  the  ions  to  the 
surface.  The  ion  current  and  energy  distribution  is  controlled  in  some  way  by  the  rf 
biasing  parameters,  the  power,  current,  frequency,  etc..  Here  we  discuss  models  of 
collisionless  rf  sheaths  in  various  regimes.  Particularly,  we  specify  a  parameter  space 
which  characterizes  rf  sheath  in  all  density  and  current  regimes,  give  the  boundaries  of 
the  various  regions,  and  give  approximate  analytic  solutions  for  the  sheath  in  each  region. 

The  standard  picture  of  an  rf  sheath  has  been  developed  by  Godyak  [1,2],  Lieberman 
[3,4],  and  others.  In  this  picture,  the  rf  frequency  is  high  compared  to  the  ion  plasma 
frequency,  so  the  ions  do  not  respond  at  the  drive  frequency.  Thus  the  theory  is 
particularly  applicable  to  capacitively  driven  plasmas  which  typically  have  rather  low 
density  (typically  less  than  1010  cm'3  and  an  rf  frequency  of  13  MHz).  In  this  theory,  as 
the  electrons  oscillate  away  from  the  bounding  surface,  a  positive  space  charge  is  left 
behind,  and  this  accelerates  the  ions  to  the  surface.  Thus  the  ions  have  a  crucial  dc 
response,  but  no  rf  response.  The  Godyak  Lieberman  (GL)  theory  is  sketched  out  in 
Section  EL  There,  in  addition  to  the  basic  theory,  the  regime  of  validity  is  specified. 

As  modem  processing  steps  put  more  demands  on  the  plasma,  usually  one  is  forced  to 
higher  density  plasmas,  with  a  more  controllable  sheath.  (In  the  capacitive  discharge,  the 
rf  creates  both  the  sheath  and  the  plasma).  Examples  of  such  higher  density  plasma 
reactors  would  be  the  electron  cyclotron  resonance  reactor  [5],  inductive  discharge  reactor 
[6],  helicon  reactor  [7],  and  large  area  plasma  processing  systems  (LAPPS)  [8].  At  higher 
density,  the  condition  for  neglect  of  ion  rf  effects  becomes  invalid.  In  fact  at  13  MHz,  the 
ion  plasma  frequency  can  be  considerably  larger  than  the  drive  frequency.  One  theory  in 
this  regime  was  developed  by  Metze,  Ernie  and  Oscam  [9]  (MEO).  This  theory 
effectively  assumes  that  the  ions  respond  to  the  rf  fields  instantaneously.  In  the  MEO 
reference,  numerical  results  only  were  presented,  so  it  is  not  clear  just  what  the  limits  of 
validity  for  the  theory  are,  and  there  were  no  simple  analytic  models.  A  more  recent  work 
[10]  modeled  the  transition  between  the  GL  and  MEO  theories,  and  results  of  numerical 
simulations  were  presented.  Again  very  little  analytic  insight  into  the  different  regimes 
was  offered.  In  Section  III  we  review  the  MEO  theory,  give  simple  analytic  models  and 
discuss  the  regimes  of  validity. 


Manuscript  approved  March  8,  1999. 
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As  a  simple  way  to  envision  the  parameter  space  for  rf  sheaths,  consider  the  single 
parameter  cox  where  co  is  the  rf  frequency  and  t  is  the  time  it  takes  the  ion  to  cross  the 
sheath.  An  ion  crosses  the  sheath  for  two  reasons,  first  because  it  enters  the  sheath  with 
some  velocity,  and  second  because  the  fields  in  the  sheath  accelerate  the  ion  across  the 
sheath.  Consider  first  the  convection  with  the  initial  ion  velocity.  The  ion  typically 
enters  the  sheath  with  the  sound  speed  cs=  [Te/M],/2  where  Te  is  the  electron  temperature 
and  M  is  the  ion  mass.  In  the  bulk  plasma,  the  electrons  carry  an  rf  current  density 
J=nevos  where  n  is  the  electron  number  density  and  vos  is  the  oscillating  electron  velocity 
in  the  bulk  plasma.  However  close  to  the  wall,  the  electrons  cannot  oscillate  in  this  way 
because  they  bump  into  the  wall.  Thus,  near  the  wall,  there  is  a  region  devoid  of 
electrons;  the  ions  crossing  this  region  in  their  self  consistent  fields  form  the  rf  sheath.  A 
first  rough  estimate  of  the  length  of  this  region  is  s0  =  vos/co.  Thus  if  the  motion  of  the 
ions  to  the  wall  with  their  input  velocity  were  the  only  consideration,  an  estimate  for  cox 
would  simply  be  vos/cs.  If  this  parameter  is  small,  the  ion  crosses  the  sheath  quickly 
because  the  sheath  is  short  and  the  ions  enter  at  high  speed.  On  the  other  hand,  if  this 
parameter  is  large,  the  ions  cross  the  sheath  slowly  compared  to  co1. 

Now  let  us  consider  the  ion  acceleration  by  the  fields  of  the  sheath.  Let  us  assume  that 
because  the  plasma  is  a  good  conductor,  the  electric  field  in  the  plasma  is  negligible.  In 
the  sheath  region,  the  residual  ions  give  rise  to  an  electric  field  of  4raies  where  s  is  the 
distance  from  the  instantaneous  edge  of  the  electrons  (is  a  function  of  time  because  the 
sheath  edge  oscillates).  Considering  the  acceleration  of  the  ions  in  this  field,  we  find  that 
the  acceleration  of  the  ion  to  the  wall  gives  a  value  of  cox  of  about  (02/(0pi2  where  tOpi  is  the 
ion  plasma  frequency.  For  small  values  of  this  parameter,  the  ions  accelerate  across  the 
sheath  very  quickly  compared  to  co'1. 

Thus  we  can  envision  a  parameter  space  characterizing  rf  sheaths  with  C0pj/<0  being  the 
horizontal  axis,  and  vos/cs  being  the  vertical  axis.  This  values  of  these  parameters 
compared  to  unity  characterizes  the  various  possible  rf  sheaths.  It  is  only  in  the  upper 
region  along  the  vertical  axis,  tOp/oxl  and  vos/cs>l  that  the  ion  remains  in  the  sheath 
period  for  a  time  long  compared  to  the  oscillation  period  (here  we  use  <  and  >  rather  than 
«  and  »  so  as  to  more  easily  show  the  regions  of  validity  on  a  two  dimensional  plot ). 
This  is  the  region  of  the  GL  sheath  model. 

Hence  the  GL  theory  is  not  only  invalid  at  high  density,  it  is  also  invalid  at  low  rf  current, 
a  fact  that  does  not  seem  to  be  generally  realized.  In  Section  IV  we  discuss  the  sort  of  rf 
sheath  found  in  this  regime.  This  sort  of  sheath  (low  density,  low  current)  has  parameters 
that  are  most  likely  not  very  interesting  for  plasma  processing  applications.  Finally  in 
Section  V,  we  enumerate  where  all  of  the  different  models  are  valid  on  the  two 
dimensional  parameter  space  we  have  specified.  Shown  there  is  this  chart  with  the 
regions  of  validity  of  the  various  sheath  models  denoted. 
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II  The  Godyak  Lieberman  [GL]  rf  Sheath  Model. 

Inside  the  bulk  plasma  there  is  a  spatially  uniform  rf  current  carried  by  the  electrons  and 
given  by 


r  t  ds 

J  =  Jn  cos  at  =  ne — 
0  dt 


(1) 


where  n  is  the  electron  density  in  the  bulk  plasma  at  the  sheath  edge  (for  simplicity  we 
assume  here  that  the  bulk  plasma  is  uniform),  and  s  is  the  position  of  the  electron. 
However  near  the  wall,  the  electron  pulls  away  and  an  ion  charge  remain  there.  In  the 
bulk  the  maximum  excursion  of  an  electron  is 


2s„ 


2JC 

neco 


(2) 


However  in  the  sheath,  the  ion  density  may  not  be  constant,  but  has  a  density  profile  n(x) 
where  s  is  the  distance  from  the  wall.  If  we  assume  that  the  plasma  is  a  good  conductor, 
and  the  electron  sheath  is  at  position  s(t),  we  have  from  Poisson’s  equation  that  the 
electric  field  in  the  sheath  is  given  by 

x 

E  =  Ane  J"  dxnt{x)  (3) 

m 


The  rf  current,  conduction  plus  displacement,  must  be  continuous  in  the  plasma  and 
sheath,  so  we  find 


ds 

n(s)e — -  J 0  cos  cot 
dt 


(4) 


where  we  have  assumed  that  the  electric  field  in  the  plasma  vanished  since  the  plasma  is  a 
good  conductor.  Thus,  if  n(s)  decreases  in  the  sheath,  as  one  goes  from  the  plasma  edge 
to  the  wall,  the  size  of  the  sheath  increases.  Now  say  the  ions  enter  the  sheath  at  smax  and 
there,  satisfy  the  Bohm  condition,  so  that  the  current  density  is  necs,  and  there  is  a 
potential  drop  of  <)>  across  the  sheath.  Now  at  t=-7c/2co,  the  electron  is  at  its  maximum 
excursion  on  the  left  (the  wall),  and  at  t=7t/2(0  it  is  at  its  maximum  excursion  on  the  right 
(the  sheath  edge).  Multiply  by  dt  and  integrate  the  right  hand  side  of  the  equation  from  - 
71/2®  to  7t/2®  and  the  left  hand  side  from  0  to  smax,  the  sheath  width.  To  integrate  the  left 
side,  express  n(s)  as  a  function  of  potential  by  using  conservation  of  energy  (assuming 
that  e<j)/Te»l)  and  assume  that  the  electric  field  in  the  sheath  is  has  the  dependence  on 
potential  as  in  a  Langmuir  Child  diode  (that  is  the  <j)(x)=<})(x/smax)4/3).  Then  we  find  that 


s 


max 


2  (2ed> 

—  - s 

3  V  T€  ° 


(5) 
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To  the  ions  entering,  the  sheath  looks  very  much  like  an  ion  diode.  However  instead  of 
vacuum,  the  ions  are  entering  a  sheath  region  where  electron  charge  oscillates  back  and 
forth,  canceling  the  ion  charge  where  it  exists.  Thus  the  average  ion  charge  is  about  half 
of  the  charge  in  a  vacuum  ion  diode.  Hence  the  simplest  theory  would  predict 

n0ecs  =  2Jlc  (6a) 

where  Jlc  is  the  Langmuir  Child  value.  For  the  current  of  an  ion  diode,  we  use  the 
expression  in  cgs  units 


7 LC  ~ 


[2e  <T 
V  M  9 ™max2 


(6b) 


on  more  familiar  units 


JLC(A/cm2)  =  10“8 


Ma  <pm (Volts) 
M  S  max2  (cm) 


(6c) 


where  MA  is  the  mass  of  the  argon  atom.  This  then  relates  the  potential  to  the  current  and 
diode  width.  A  schematic  of  the  GL  sheath  model  is  shown  in  Fig.  (1). 

The  more  precise  theory  of  Lieberman  [3,4]  predicts  a  factor  of  1.8  instead  of  2  in  Eq. 
(6a).  At  his  point  it  might  appear  that  there  is  a  dc  (ion)  current  to  the  sheath.  However 
this  is  not  necessarily  so.  Recall  that  the  electron  sheath  oscillates  in  time  between  the 
maximum  position  and  the  substrate.  When  the  electron  sheath  is  in  contact  with  the 
substrate,  an  electron  current  can  flow  to  the  substrate  which  is  essentially  unlimited. 

Thus  when  the  electron  sheath  hits  the  surface,  the  dc  current  can  be  canceled  or  can  take 
on  whatever  value  is  imposed  by  the  external  circuit. 


In  any  case,  solving  the  modified  Langmuir  Child  diode  law,  Eq.  (6),  for  the  voltage,  and 
using  the  fact  that  the  sheath  width  now  depends  on  voltage  according  to  Eq.  (5),  we  find 


<p(Volts)  = 


0A5Jg\mA!  cm2) 


T(ev)\ 


(7) 


where  f  is  the  rf  frequency.  This  gives  a  simple,  approximate  law  for  the  dc  sheath 
voltage  in  terms  of  the  plasma  and  rf  parameters.  Notice  that  the  calculation  of  the  sheath 
width  smax  depends  on  the  functional  form  of  <}>(x)  in  the  sheath  region.  The  assumption 
of  a  relation  appropriate  for  a  space  charge  limited  diode  gives  the  coefficients  in  Eqs.(5 
and  7).  Other  functional  forms  for  <[>(x)  would  give  different  coefficients  in  Eq.(5).  For 
instance  the  assumption  of  a  uniform  electric  field  would  give  an  smax  larger  by  a  factor  of 
1.5.  However  this  factor,  when  raised  to  the  fourth  power  would  give  a  voltage  larger  by 
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period  2rc/co. 


a  factor  of  5  in  Eq.(7).  Thus  the  voltage  is  difficult  to  predict  accurately  from  simple 
considerations;  however  the  sheath  width  should  be  predicted  reasonably  accurately.  It 
might  be  best  to  regard  Eq(7)  as  a  scaling  law,  with  a  coefficient  to  be  determined 
empirically. 


There  are  two  conditions  under  which  the  GL  model  is  valid.  First  of  all  the  ions  are 
assumed  to  not  react  to  the  rf  fields,  but  only  the  dc  fields.  Clearly  this  implies  C0pj<C0,  or 
the  ions  can  respond  to  the  rf  as  well  as  dc  fields.  (In  the  bulk  plasma  the  ion  dielectric 
constant  C0pj2/(02  is  simply  added  to  the  much  larger  electron  dielectric  constant  (Ope2/co2,  so 
the  ions  do  not  play  a  significant  rf  role  in  the  bulk,  at  least  for  an  unmagnetized  plasma.) 
The  second  condition  is  that  the  ion  transit  time  across  the  sheath  is  very  long  compared 
to  2tc/co,  so  that  the  ions  only  react  to  the  average  fields  in  crossing  the  sheath.  The  sheath 
width  is  about  [2e(J)/Te]  1/2s0,  and  the  ion  velocity  across  the  sheath  is  about  [2e<j)/M]1/2, 
where  M  is  the  ion  mass.  Hence  the  time  for  the  ions  to  cross  the  sheath  is  about  So/cs.  In 
order  for  this  to  be  large  compared  to  the  oscillation  time,  we  must  have  vos/cs>7t,  where 
vos  is  the  oscillating  velocity, 


If  this  condition  is  violated,  the  ions  react  to  the  time  oscillation  of  the  sheath  and  there 
will  be  a  spread  in  ion  energies  impinging  on  the  substrate.  Notice  that  the  GL  model 
predicts  monoenergetic  ions  impinging  on  the  substrate.  The  monoenergetic  nature  of  the 
ion  flux  can  be  disrupted  by  either  ion  rf  effects,  or  ion  transit  time  effects.  Thus  the  two 
conditions  for  the  validity  of  the  GL  model  of  the  sheath  are 


— <  1 
co 


(9a) 


v 

-^>1 

Cs 


(9b) 


If  we  envision  a  parameter  space  where  the  horizontal  axis  is  cfyj/co  and  the  vertical  axis 
is  vos/cs,  then  the  GL  model  of  the  rf  sheath  is  valid  along  the  vertical  axis  above  vos/cs>l. 
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Ill  The  Metze,  Ernie,  Oskam  rf  Sheath  Model 

The  MEO  model  assumes  that  the  ions  cross  the  sheath  as  if  the  sheath  were  constant  in 
time.  This  means  that  in  all  cases,  cot<l.  However  the  ion  plasma  frequency  must  be 
large  compared  to  co  so  that  ion  fluid  responds  to  the  instantaneous  fields  of  the  plasma. 
The  sheath  current  is  assumed  to  consist  of  both  displacement  current  and  conduction 
current.  Thus  there  are  two  possible  regimes  for  the  MEO  sheath  one  where  conduction 
current  dominates  and  one  where  displacement  current  dominates.  We  will  devise  simple 
analytic  models  for  each.  The  Poisson’s  equation  for  the  sheath  at  any  time  is  given  by 

44  =  — 4t n?(n(.  -  ne )  (10) 

ax 

For  the  electron  and  ion  densities  we  assume  the  normal  formulae  for  Langmuir  sheath 
theory. 


f 

ne  =  «ex  p 


(Ha) 


n 


(lib) 


Here  we  have  assumed  that  where  each  density  is  equal  to  n,  the  ion  velocity  into  the 
sheath  is  cs  and  the  potential  at  this  point  vanishes.  Going  further  into  the  sheath,  the 
potential  gets  more  and  more  negative,  accelerating  ions  into  the  wall  and  repelling 
electrons.  Integrating  Eq.(10)  once  by  quadrature,  we  find  the  electric  field  in  terms  of 
the  potential  to  be 


E  =  -JSmT( 


2e(f> 


+  exp 


-11/2 


-2 


V  -W 


(12) 


where  the  sign  of  E  is  chosen  such  that  ions  are  accelerated  into  the  wall  on  the  left,  as  in 
Fig.  (1).  Thinking  of  a  potential  drop  across  the  sheath  which  is  considerably  larger  than 
the  electron  temperature,  we  can  adopt  the  following  approximation  for  E  in  terms  of  <|) 


E  =  -yfSmf( 


E  =  0 


<(k0 

(13a) 

(j»0 

(13b) 
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From  Poisson’s  equation,  we  can  determine  that  the  surface  charge  density  o  on  the 
substrate  is  simply  equal  to  E/4jt,  so  one  can  derive  a  charge  in  terms  of  sheath  potential. 
Defining  the  capacitance  per  unit  area  as  C’=do/d<}>,  we  find  that 


C  =  0.25 


\nTe 

2e 

iln 

_  <i>X. 

<jx0 


(14a) 


C’=0 


(t»0 


(14b) 


Thus  one  equation  relating  the  sheath  voltage  and  current  is 


C  —  +  7=7 
dt  c  * 


(15) 


Here  Jc  is  the  conduction  current,  and  Jrf  is  the  rf  current  J0coscot.  The  first  term  is  of 
course  the  displacement  current.  There  are  two  possible  limits  to  Eq.  (15),  either  the 
conduction  current  dominates  or  the  displacement  current  dominates. 


We  consider  first  the  former.  The  conduction  current  as  a  function  of  substrate  voltage  is 
as  given  for  a  Langmuir  probe, 


Jc  =  ne[Te  /  2mn]exp(e(j)  /Te)~  necs 


(16) 


where  m  is  the  electron  mass,  and  J  is  plotted  as  a  function  of  <J)  in  Figure  (2).  For 
positive  potential,  there  is  the  large  electron  current,  while  for  negative  potential  the 
current  density  asymptotes  to  a  much  smaller  negative  value,  the  ion  saturation  current 
necs.  The  actual  dependence  of  J  on  <J)  for  positive  <J>  is  not  important  here,  only  the  fact 
that  it  is  large  and  rapidly  increasing.  Actually  the  current  density  asymptotes  for 
negative  values  to  necs,  and  the  voltage  is  determined  by  the  input  power.  For  the  other 
half  cycle  the  current  is  electron  current,  and  the  sheath  voltage  is  small,  as  is  the  input 
power.  The  ions  impinging  on  the  substrate  have  a  spread  in  energy  from  about  zero  up 
to  <j>max?  determined  by  the  input  power.  The  energy  of  a  particular  ion  depends  on  what 
time  it  enters  the  sheath. 


We  now  estimate  the  sheath  potential  in  terms  of  the  input  power.  To  do  so,  we 
approximate  the  graph  of  Fig.  (2)  as  that  of  a  diode  with  back  leakage  current  density 
equal  to  -necs.  For  such  a  diode,  connected  to  an  external  circuit,  with  negative  voltage, 
the  voltage  is  then  specified  as  a  function  of  time.  In  the  negative  half  cycle,  we  assume 
that  <J)(t)=-<|)sin  cot  for  0<t<rc/co.  During  this  half  cycle,  the  current  density  is  -necs.  During 
the  positive  half  cycle,  both  the  current  and  voltage  will  be  small,  the  main  constraint 
being  that  the  total  dc  current  in  the  entire  cycle  is  zero  (or  else  is  as  specified  by  the 
external  circuit).  In  terms  of  the  power  per  unit  area  going  into  the  sheath,  Pcx,  one  can 
determined  the  sheath  voltage.  The  result  is 
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(17) 


Pcx  =  —<pnecs 
71 

where  we  have  assumed  that  there  is  no  power  dissipated  on  the  positive  half  cycle. 

There  are  several  conditions  which  must  be  satisfied  in  order  to  have  a  conduction  current 
dominated  rf  sheath.  First  of  all  the  maximum  rf  current,  J0,  can  only  be  equal  to  necs,  or 

-^-<1  (18a) 

c. 

Secondly,  the  conduction  current  must  dominate  the  displacement  current.  This  gives  rise 
to  a  maximum  frequency.  Denoting  the  maximum  potential  by  <(),  and  comparing  the 
conduction  current  to  the  displacement  current,  we  find  that  conduction  current  dominates 
if 


2.5 


V  r  ) 


1/4 


(0  ) 


\e<Pj 


>1 


(18b) 


Unless  the  potential  is  very  large,  this  is  nearly  equivalent  to  the  to  the  condition  for  the 
ions  to  cross  the  sheath  quickly. 


to 


PL>[ 

a ) 


(18c) 


These  define  the  region  for  an  MEO  sheath  with  conduction  current  dominant  (we  call 
this  an  MEOCC  sheath).  In  the  parameter  space  defined  at  the  end  of  Section  II,  the 
MEOCC  sheath  exists  along  the  horizontal  axis  for  values  of  tOpj/co  Greater  than  unity. 


We  now  consider  The  MEO  sheath  where  displacement  current  dominates  We  call  this 
the  MEOco  sheath.  The  equation  for  this  is 


2e  ' 

J  2n 

.  4>X_ 

1/4 


dl 

dt 


=  -J„  sin  cot 


(19) 


Of  course  the  displacement  current  can  only  dominate  for  J<0,  so  we  consider  the 
solution  to  Eq.(19)  for  times  0<t<7t/co.  Let  us  say  that  at  t=0,  <j)=0.  Then  we  find  that  the 
potential  as  a  function  of  time  is 


-m= 


32tc2J0 4 1"  cos  cot  - 1 
H27>ft)4  L  2 


COS  (01  - 1  ~|4 
- 


0<t<7t/O) 


(20a) 


or  an  expression  for  the  maximum  voltage  in  more  familiar  units  is 
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(20b) 


( p(volts )  =  2.4 
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The  scaling  of  the  Voltage  with  parameters  is  the  same  as  in  the  GL  sheath,  but  the 
coefficient  is  larger  by  about  a  factor  of  5.  However  the  ions  do  not  strike  the  wall  with  a 
mono-energetic  energy  distribution  any  longer,  depending  on  the  time  they  enter  the 
sheath,  their  energy  varies  from  zero  to  the  maximum.  The  average  value  of  the  ion 
energy  is  three  eighths  the  maximum  value.  Furthermore,  this  ion  acceleration  exists  only 
for  half  of  the  rf  cycle  (we  will  discuss  the  other  half  cycle  shortly).  Hence  a  time  average 
over  the  entire  cycle  reduces  the  effective  potential  by  an  additional  factor  of  2,  so  that 
when  all  the  averaging  is  done,  potential  is  reasonably  comparable  to  that  in  the  GL 
model. 


Let  us  now  consider  the  dynamics  of  the  formation  of  the  rf  sheath  in  the  regime  where 
the  displacement  current  dominates.  At  t=0,  the  plasma  is  in  contact  with  the  wall  when  a 
negative  voltage  is  applied.  Ions  are  accelerated  into  the  wall  and  the  electrons  move 
away.  At  each  instant,  the  sheath  region  appears  to  be  a  Langmuir  Child  ion  diode  as 
specified  by  Eq.  (6).  However  now  the  current  is  the  displacement  current, 

ds 

J  —  ne —  (21) 

dt 


For  a  voltage  pulse  specified  as  a  function  of  time,  we  can  multiply  the  Langmuir  Child 
relation,  Eq.(6)  by  s2  and  integrate  each  side  in  time  from  zero  to  the  voltage  maximum. 
The  left  hand  side  simply  gives  s3/3.  This  is  the  basic  idea  used  in  calculating  the  ion 
dose  to  a  substrate  in  plasma  immersed  ion  implantation  [1 1,12].  Assuming  a  voltage 
pulse  given  by  Eq.(20),  we  find  that 

2.6*10  ) 

»_(cm)  =  — 7 - v,v  7  ^  y  (22) 

(bO  (life) 

Assuming  that  the  ions  accelerates  through  the  sheath  through  the  potential  <J>,  we  find  that 
arc  is  given  by 


cor  ~0.4 

\ 


CO 


\  2/3 


(23) 


confirming  our  original  condition  for  the  validity  of  the  rf  sheath  in  this  region.  Thus  the 
MEO  sheath  with  conduction  current  dominant  occurs  for 
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that  is  the  upper  region  in  the  parameter  space  we  have  specified.  Once  the  voltage 
reaches  its  maximum,  the  current  changes  sign,  and  the  rf  sheath  quickly  collapses.  For 
positive  current,  both  conduction  and  displacement  current  may  be  important,  but  in  any 
case  the  voltage  is  low,  and  the  total  current  in  this  half  cycle  simply  cancels  the  total 
current  in  the  first  half  cycle. 

There  is  one  additional  problem  regarding  the  MEOco  sheath.  Every  cycle  the  sheath 
reaches  into  the  plasma  and  pulls  nsmax  ions  per  unit  area  into  the  wall.  Thus  the  average 
flux  of  ions  is  nsmaxco/27t.  However  the  maximum  flux  which  can  be  supplied  by  the 
plasma  external  to  the  sheath  is  ncs.  Here  recall  that  n  is  the  value  at  the  edge  of  the 
sheath.  Thus  the  MEOco  sheath  can  only  function  as  described  if 

C  ( I ) 

-****— <1  (25a) 

2  KC, 


Inserting  the  parameters  for  smas  which  we  have  derived,  we  find  that  in  the  cOp/co,  vos/cs 
plane,  this  condition  reduces  to 


\Cs  J 


® eL 
CO 


,4/3 
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(25b) 


Thus  in  our  parameter  space,  in  the  upper  most  region,  the  plasma  is  not  able  to  supply 
the  necessary  flux  demanded  by  the  rf  sheath.  In  this  upper  region,  for  the  rf  sheath  to 
exist,  it  appears  the  rf  must  modify  the  ambient  plasma.  It  may  do  this  in  at  least  three 
ways.  First  it  can  supply  additional  electron  heating  so  as  to  appropriately  increase  cs; 
secondly,  it  may  drain  the  plasma  until  the  density  is  low  enough  that  Eq.(25b)  is 
satisfied;  and  thirdly,  it  may  produce  additional  ionization  in  the  sheath  so  that  the  ion 
drain  is  equal  to  the  ion  input  (flowing  in  at  the  sound  speed)  plus  the  ions  produced  in 
the  sheath  region  by  the  rf  fields.  In  other  words,  the  rf  sheath  not  only  accelerates  the 
existing  plasma  to  the  wall,  it  also  must  modify  the  plasma  in  the  sheath  region  and/or  in 
the  bulk.  The  region  of  our  parameter  space  where  the  sheath  does  not  modify  the  plasma 
is  rather  limited,  and  considering  the  approximations  made,  the  MEOco  sheath  without 
plasma  modification  may  exist  only  as  an  idealization. 


IV  The  Modified  DC  (MDC)  Sheath 

Let  us  now  consider  the  rf  sheath  in  the  case  of  low  current  and  low  density.  By  low 
density,  we  mean  that  the  ions  are  unable  to  respond  the  rf  fields,  and  only  respond  to  the 
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dc  fields.  This  means  of  course  that  (Opi/oxl.  Also  we  assume  that  vos/cs<l,  so  as  to 
distinguish  it  from  the  low  density  Godyak/Lieberman  model  of  an  rf  sheath.  At  low 
density  and  low  current,  a  matrix  sheath  model  [13]  is  a  reasonable  one;  that  is  a  model  in 
which  the  electrons  oscillate  away  from  the  wall,  but  the  ion  charge  density  is 
unperturbed.  The  width  of  the  sheath  is  now  2so=2Jo/ne(0,  the  maximum  potential  drop 
between  the  sheath  edge  and  wall  is 


<j)  =  limes2  =  Sime 


(26) 


2 

and  the  maximum  ion  energy  gain  compared  to  the  entrance  energy,  0.5Mcs  is  given  by 


A £„ 


\Mc- 


4 O)  2v  2 

pi  os 
2 

CD  Cs 


(27) 


Thus,  in  the  regime  in  which  we  have  defined  the  sheath,  the  ion  energy  change  is 
relatively  small.  This  justifies  the  assumption  that  ion  density  is  unperturbed.  The  time 
for  an  ion  to  cross  the  sheath  is  roughly  s,Jc%.  Comparing  this  to  to’1,  we  see  that  the 
crossing  time  is  smaller  than  to 1  by  a  factor  of  vos/cs. 


However  there  is  normally  a  dc  sheath  if  the  wall  is  at  the  floating  potential  (J=0).  Thus 
the  dc  sheath  has  an  additional  negative  potential  drop  of  0.5Teln[M/27tm].  Hence  from 
the  dc  sheath  alone  the  ions  strike  the  wall  with  energy  0.5Te{  1  +  ln[M/2jrm]}.  The 
effect  of  the  rf  is  then  only  to  slightly  modify  the  dc  sheath,  by  spreading  the  energy  of 
ions  striking  the  wall  over  an  energy  range 
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In  the  parameter  space  we  have  defined,  this  sheath  is  characterized  by  the  region  near  the 
origin.  Since  it  only  modifies  the  dc  sheath  slightly,  and  in  doing  so,  only  spreads  out  the 
nearly  uniform  energy  distribution,  one  would  not  normally  use  an  rf  sheath  in  this 
regime  for  plasma  processing. 


V  Parameter  Space  for  Collisional  rf  Sheaths 

We  have  seen  that  there  are  a  number  of  different  possible  regimes  for  the  behavior  of  rf 
sheaths  in  collisionless  plasmas.  It  is  possible  to  plot  these  regions  on  a  simple  two 
dimensional  plot  where  the  horizontal  axis  is  (Opj/co,  and  the  vertical  axis  is  vos/cs.  Figure 
3  shows  such  a  plot.  We  have  found  simple  analytic  models  of  the  sheath  in  each  of  the 
regions.  Also  shown  in  the  MEOco  region  is  the  curve,  above  which  the  sheath  cannot 
exist  unless  the  rf  which  forms  the  sheath  also  increases  the  electron  density  and/or 
electron  temperature  at  the  sheath  edge,  so  that  the  exterior  plasma  can  supply  the 
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necessary  particle  flux.  The  basis  of  these  models  are  that  the  ions  density  is  either  so 
low,  that  the  ions  cannot  respond  to  the  high  frequency  fields;  or  is  so  high  that  they 
respond  to  the  rf  fields  as  if  they  were  dc.  The  only  model  which  predicts  a 
monoenergetic  flux  of  ions  impinges  on  the  substrate  is  the  low  density,  high  current 
limit,  the  GL  sheath.  The  ions  pick  up  a  distribution  in  energy  as  one  goes  to  either  high 
density  or  low  current. 
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